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Abstract 

Thermodynamical fluctuations of temperature in mirrors of gravitational 
wave antennae may be transformed into additional noise not only through 
thermal expansion coefficient ||] but also through temperature dependence of 
refraction index. The intensity of this noise is comparable with other known 
noises and must be taken into account in future steps of the antennas. 

I. INTRODUCTION 

We have shown in our previous article that thermodynamical fluctuations of tempera- 
ture in mirrors (test masses) of LIGO-type gravitational wave antenna are transformed 
due to the thermal expansion coefficient a = {l/l){dl/dT) into additional (thermoelastic) 
noise which may be a serious "barrier" limiting sensitivity. This noise is caused in fact 
by random fluctuations of the coordinate averaged over the mirror's surface. The spectral 
density of this random coordinate displacement may be presented for inflnite test mass in 
the following form 



^ This result was refined for the case of finite test masses by Yu. T. Liu and K. S. Thorne 
However the difference from our calculation is only several percents for the planned sizes of test 
masses, and hence we use here much more compact expression (|l|). 
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Here k is the Boltzmann constant, T is temperature, a is Poison ratio. A* is thermal con- 
ductivity, p is density and C is specific heat capacity, ro is the radius of the spot of laser 
beam over which the averaging of the fluctuations is performed. This noise is of nonlinear 
origin as the nonzero value of a is due to the anharmonisity of the lattice. 

The goal of this article is to present the results of the analysis of another additional (and 
also of nonlinear origin) effect which may be comparable with other known noise mechanisms 
limiting the sensitivity. Qualitatively this effect is may be understood in the following way. 
The laser beam "extracts" the information not only about the change of the length of 
the antenna produced by gravitational wave but also about the fluctuations of position of 
mirrors' surfaces averaged over the beam spot. These fluctuations lead to phase noise in the 
reflected optical fleld. However the phase noise may be produced by another effect. High 
reflectivity of the mirrors is provided by multilayer coatings which consist of alternating 
sequences of quarter- wavelength dielectric layers having refraction indices rii and 77,2. The 
most frequently used pairs of layers are Ti02 — Si02 and Ta20z — Si02- While reflecting * 
the optical wave "penetrates" in the coating on certain depth. This depth is of the order 
of the optical thickness of one pair of layers I < Ifi. If the values of ni and 77-2 depend 
on temperature T (thermo-refractive factor j3 = dn/dT is nonzero) then thermodynamical 
fluctuations of temperature lead to fluctuations of optical thickness of these layers and hence 
to the phase noise in the reflected wave. Though the thickness / of the working layer is small, 
the coefficient j3 is usually significantly larger than a (both have the same dimensions). For 
fused sihca {Si02) a = 5 x 10"''' and /3 = 1.45 x 10~^ (i.e. 30 times larger than a). 
This phase noise may be evidently easily recalculated in terms of equivalent fluctuations of 
the surface and consequently compared with the spectral sensitivity of the antenna. 

We have analyzed also the photo-thermal refractive shot noise: due to random absorption 
of optical photons, the random fluctuations of temperature in the surface layer of the mirror 
take place, producing fluctuations of refractive indices of the coating and therefore phase 



2 



fluctuations of reflected light wave (this effect is similar to photo-thermal shot noise, analyzed 
in However, this effect is numerically much smaller than thermo- refractive noise — that 
is why we do not present here the detailed analysis of it. 



II. THERMO-REFRACTIVE NOISE 

The theory of reflection of light from multilayer dielectrical coating is well known (see 
for example [^). Using traditional approach we may recalculate the phase shift 50 into 
equivalent displacement 5x of mirror (see Appendix ^): 

^ ^ -^*A/9efr, (2) 

~ Ainl - nl) ■ ^"^^ 
Here u is the fluctuation of averaged temperature, jSi = drii/dT, (32 = dn2/dT . It is 
important to note, that effective coating thickness is much smaller than the characteristic 
length of diffusive heat transfer: / ^ a/ ^Juj [a is temperature conductivity, uj is the frequency 
of observation which is of order ~ lOOi/z for laser gravitational wave antenna). Therefore 
we may consider in our calculations that fluctuations of temperature are correlated in the 
layers. 

To calculate thermodynamical fluctuations of temperature M(r, t) in the surface layers 
we use Langevin approach and introduce fluctuational thermal sources F{f,t) added to the 
right part of the equation of thermal conductivity: 

^-a^Au = Fif,t), a' = ^. (4) 

This approach was described and verifled in (see all the details over there). As in 
we replace the mirror by half-space: — oo<x<oo, — oo<y<oo, < z < oc with 

the boundary condition of thermo-isolation on surface z = 0. We may now calculate the 

spectrum of temperature fluctuations: 

— * 
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{F{k, u;)F*{k,, uj,)) = (27r)^|A;p 6(00 - 00,) x (7) 

X 6{k^ - k^i)6{ky - kyi) x 
X [5{k-,-k-,i) + 5{k^ + k-,i)]. 

The thermodynamical fluctuations of temperature u averaged over tlie volume V = Tir^l 
may be presented in tlie following form: 



u 

^'0 



\- r dxdy r dz u{f,t) e~("'+^')/'^o e-^/' 



TTTnl J -00 Jo 

dkduj F{k,Uj)e''^' ^_(;t2+fc2),2/4 1 



(27r)4 a^\k\'i + iuj l-ik,V ^ ' 

From this expression and from (|^) we find immediately the spectral density S'„(co') of 
fluctuations of the averaged temperature: 

^ , , 2kT2 a* 27r k^dkx. dk. 



(pC)2 Jo (27r)2 7-00 2n 

kl + fci ,-feir2/2_(l + 1) 
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a4(A;2 + fc2)2 + a;2 1 + A;2/2 
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Here k\ = k^. + ky. The first term 2 appears because as in we use "one-sided" spectral 
density, defined only for positive frequencies, which is connected with the correlation function 
{u{t)u{t + r)) by the formula Su{uj) = 2J^^dT {u{t)u{t + r)) cos(cjr). The term (1 + 1) 
appears due to two (5-functions in square brackets in (0). For the frequency of observation 
^ 27r X 100 s~^ characteristic length a/^/cJ ~ 50 /i (we used for the estimates constants 
for fused silica), so that / <^ a/^/cJ ^ tq. Taking into account that A;^ ~ 1/ro ^ ^/oj/a we 
may consider the first denominator as constant while integrating over In the same way 
kz — l/l ^ \/Toja and while integrating over k^ we may consider the second denominator as 
unity. It is interesting that in this model the spectral density ^^(ti;) does not depend on /. 

This spectral density may be recalculated to the spectral density of equivalent fluctua- 
tions of surface displacement to compare it with other known sources of noise: 
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FIGURES 

Noises in gravitational wave antennae 
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FIG. 1. Comparison of SQL-limited sensitivity with different sources of noise in gravitational 

f, (Hz) 

wave antennae: thermo-refractive, Brownian (dominating in fused silica mirrors) and thermo-elastic 
(dominating in sapphire mirrors). 



III. NUMERICAL ESTIMATES 



For the numerical estimates we assumed that the multilayer coating consists from alter- 
nating pairs of layers: Ti02 {ni = 2.2) and Si02 (^^2 = 1-45), or Ta20^ {rii = 2.2) and 5*^02 
(^2 = 1.45). The values of f3 for Ti02 and for Ta20r, were found in 

We want now to compare the thermo-refractive fluctuations with thermoelastic noise (|I]) 
and noise associated with the mirrors' material losses described in the model of structural 
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damping (we denote it as Brownian motion of the surface). In this model the angle of 
losses (p does not depend on frequency and for its spectral density the following formula is 
valid for infinite test mass [fZlpI 



o 



where E is Young modulus, and a is Poison ratio. 

The spectrsal sensitivity of gravitational wave antenna to the perturbation of metric 
h{uj) may be recalculated from noise spectral density of displacement x using the following 
formula: 



h{u) = , (12) 

where we used the fact that antenna has two arms (with length L) with two mirrors the 
fluctuations on which are averaged over the radii roi and ro2. 

The LIGO-II antenna will approach the level of SQL, so we also compare the noise limited 
sensitivity to this limit in spectral form M: 



V mw^L^ 

For the calculations we used the set of parameters given in Appendix ^ (the same as in 
i) pins i 

roi = 3.6/ cm, = 4.6/-\/2 cm. 



ni = 2.2, /?2 = 4 ■ 10-5 (TzOa), 
ni = 2.2, /92 = 6 ■ 10-5 K'^ (TaaOs), 
ns = 1.45, A = 1.5 ■ 10-5 K-^ {S1O2), 



We used figures from ^ for ion plating method only, for other methods of deposition the 
value of (3 may be two times larger. In figure 1 we plot all previously known noises 
together with the new one. We see that thermorefractive noise limitation is close to SQL 
for the frequences near 200 Hz. 



CONCLUSION 



Summing up, we may say that thermo-refractive effect is not small and it must be seri- 
ously considered in interferometric gravitational antennae (projects LIGO-II and especially 
LIGO-III, where overcomming of the SQL is planned). It is also important that this ef- 
fect depends slower on the radii of the beam-spots than thermo-elastic noise and thus may 
become dominating for larger tq planned in LIGO-II and LIGO-III. 
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APPENDIX A: 

In this appendix we give the calculation of coefficient of reflection of light wave from 
multilayer coatings consisting of infinite sequences of pairs of quarter-wavelength dielectrical 
layers rii and n2. 

Let the refraction index of odd layers fluctuates on Arii and the refraction index of even 
layers on An2. One may reformulate this problem into the problem for distributed long line 
ip. The equivalent impedance Z of this system of layers may be deduced using the following 
statement: the addition of two layers does not change the value of Z. 

Voltage V2 and current I2 at the end of second layer may be found from input voltage 
Vo and current Jq using transformation matrix M ( ||^, formula (3.9.27)): 











^ Mu 


Mi2^ 


= M X 


, M = 










^M2l 


M22 J 



Mil = cos 01 cos 02 sin 0i sin 02, 

n2 



8 



/sin cos 02 sin cos 0i \ 

Mi2 = -I \ , 

\ rii n2 J 

M21 = —i {n2 sin 02 cos 0i + rii sin 0i cos 02) 

77-2 

M22 = cos 01 COS 02 sin 0i sin 02, 



M 



ni 
n2 



(^£1 A- !i 
\ni 1 



."Hi 
ni 



J 



^ i {n2(pi + ni(p2) 

Here we take into account that for quarter- wavelength layers 0i = 7r/2 + ipi, 02 = 
7t/2 + ip2 and therefore one may use approximation sin0i ~ 1, sin 02 ~ 1, cos0i ~ 

-ipl, COS02~-V52- 

Now we put that Jq = YVq and I2 = YV2 {Y = 1/Z is generalized conductivity of the 
sequence of layers) and obtain two equations: 



n2 \ni n2/ 



YV2 = Vo[i (n2V?i + 111(^2) - Y 



n2 
ni 

Solving these equations we find conductivity Y and reflection coefficient K: 

nin2 



(Al) 
(A2) 



F ~ — i 



(n2V5i + niV92) 



K 



Y - 1 

Y + 1 



nin2 

-1 - 2i 2 (^aV^i + '^iV^a) 

nf — ni 



(A3) 
(A4) 



From this point it is easy to obtain (HH), assuming 



^1 
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APPENDIX B: PARAMETERS 



u; = 27r X 100 s-\ T = 300 K, 

m = = 3 X 10^ g, A = 1.06 /i, L = 4 x 10^ cm; ; 

Fused silica: 

a = 5.5 X 10-^ K"\ A* = 1.4 x 10^ — 

cm s K 
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p = 2.2 C = 6.7x 10^ 

cm-^ g K 

E = 7.2 X 10^^— ^, (7 = 0.17, = 5x10"^; 



Sapphire: 



a = 5.0 X 10"^ K-\ X* = 4.0 x 10*^ 
p = 4.0 C = 7.9x 10^ 



erg 



cm s K' 



E = Ax 10^2 ^ = 029, = 3 X 10"^ 



erg 



cm"' 
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